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Resume — L'experience acquise avec de nombreuses applica- 
tions reussies permet de revisiter la commande sans modele. 
On remplace la derivation numerique de signaux bruites 
par une identification parametrique en temps reel, beau- 
coup plus simple. On explique la curieuse universalite des 
PID usuels et pourquoi un modele ultra-local d'ordre 1 suf- 
fit presque toujours. On montre que, meme avec un modele 
restreint, I'usage d'un correcteur PI intelligent reste avanta- 
geux. Deux exemples, un pendule inverse et un echangeur de 
vapeur, semblent conflrmer I'avantage de la commande sans- 
modele par rapport a d'autres approches, comme la com- 
mande par regimes glissants et I'utilisation de techniques de 
dimension inflnie, qu'il s'agisse de retards ou d'equations aux 
derivees partielles. Plusieurs simulations numeriques illus- 
trent notre propos. 

Abstract — The experience gained with numerous success- 
ful applications permits to revisit some points of model-free 
control. The numerical differentiation of noisy signals may 
be replaced by a real time parameter identification which is 
much simpler. The strange ubiquity of classic PIDs is explai- 
ned as well as the almost universal utilization of ultra-local 
models of order 1. We show that even with a partially known 
model the utilization of an intelligent PI controller remains 
profitable. Two examples, an inverted pendulum and a heat 
exchanger, seem to confirm the superiority of model-free 
control with respect to sliding modes control and techniques 
stemming from infinite-dimensional systems, such as delays 
or partial differential equations. Our paper is illustrated by 
several computer simulations. 



Mots-cles — Commande sans modele, PI intelligent, PID, 
identification, pendule inverse, echangeur de chaleur, re- 
gimes glissants, retards, equations aux derivees partielles. 
Keywords — Model-free control, intelligent PI, PID, identi- 
fication, inverted pendulum, heat exchanger, sliding modes, 
delays, partial differential equations. 



I. Introduction 

« Rien de plus utile qu'une bonne tlieorie0», voila qui 
resume les ambitions de la commande sans modele, due 
a deux des auteurs ([17], flSJ). Cette vision, radicalement 
nouvelle de I'automatiquelj, combine, rappelons-le, les deux 
avantages suivants : 

1. Inanite, comme son nom I'indique, d'une modelisation 
mathematique precise, qu'elle soit basee sur des lois phy- 
siques ou des procedures d'identification, comme I'exige 
trop souvent la theorie du controle « moderne ». Frotte- 
ments, hysteresis, effets thermiques, vieillissement, disper- 
sion des caracteristiques due a la fabrication en serie, . . . , 
sont, on ne le salt que trop, rebelles a toute description 
simple et fiable par equations difFerentielles. 

2. Facilite du reglage du correcteur correspondant, dit PI 
intelligent, ou iPI. D'ou une rupture nette avec les PID clas- 
siques dont le reglage, malaise et penible (voir, par exemple, 
[5j, [38J), releve trop souvent de « recettes de cuisine ». 
De nombreuses applications (voir [1], [3], [9], [26], [27], [32] . 
ES], [37], [47], [11], [49], p]) out deja ete reussies. Get 
article tire parti de l'experience acquise pour ameliorer la 
presentation. 

Voici les modifications essentielles : 
1. On remplace, comme en |20j, la derivation de signaux 
bruites ([21], [35]) par identification parametrique li- 
neaire, developpee en [53], [5S], conceptuellement plus 

1. Traduite de I'anglais : « nothing is as practical as a good 
theory », cette phrase est due au psychologue Kurt Lewin. Les au- 
teurs remercient chaleureusement Roger Tauleigne (ECS, ENSEA, 
Cergy-Pontoise) de la leur avoir suggeree comme titre. 

2. Voir, cependant, [8]. [28]. 



simple, qui a suscite par ailleurs quelques applications ([2], 

0. m, m)- 

2. On explique 

- pourquoi un modele ultra-local d'ordre 1 sufSt presque 
toujous en pratique ; 

- I'etrange ubiquite des PID classiques, malgre leurs re- 
glages fastidieux (voir [4J). 

3. On niontre, suivant [5D] et contrairement a [T7], |18j . 
qu'avec un modele restreint, c'est-a-dire partiellement 
connu, on peut conserver un correcteur iPI. 

Apres des rappels au § |lTl on traite ces questions aux § 
inil HYI El et IVD Les illustrations numeriques du § IVIH 
c'est-a-dire un pendule inverse et un echangeur de chaleur, 
mettent en evidence la probable superiorite de notre ap- 
proche par rapport aux modes glissants (voir, aussi, |20) . 
|41| ) et questionne I'utilite de techniques issues de la di- 
mension infinie, comme les retards (voir, aussi, [17], |18) ) 
et les equations aux derivees partielles (voir, aussi, [33] et 
[T]). La conclusion du § IVIIII questionne le futur de la re- 
cherche universitaire en automatique et, plus generalement, 
la place preponderante de la modelisation mathematique en 
sciences appliquees. 

II. Rappels sur le sans-modele 

A. Le modele ultra-local 

On se restreint, pour simplifier les notations, a un sys- 
teme monovariable, d'entree u et de sortie y. Dans I'igno- 
rance d'un modele mathematique global, on introduit le 
modele « phenomenologique » , dit ultra-local, valable sur 
un court laps de temps. 



+ (1) 



ou 

- I'ordre de derivation v, en general 1, choisi par I'ope- 
rateur, est etranger a I'ordre de derivation maximum 
de y, inconnu, dans le systeme ; 

- le parametre constant a, fixe par I'operateur afin que 
les valeurs numeriques de au et y^'^'^ aient meme ordre 
de grandeur, n'a pas a priori de valeur precise ; 

- F, qui contient toutes les informations « structu- 
relles » , depend de toutes les autres variables du sys- 
teme, y compris des perturbations, et de leurs derivees. 

L'estimation en temps reel de la valeur numerique de F, 
traitee au § IIII-Al permet de reactualiser ([1]) a chaque ins- 
tant. 

B. Correcteurs PI intelligents 

On obtient le comportement desire avec, si = 1 en P]). 
le correcteur proportionnel-integral intelligent, ou iPI, 



F -y* + Kpe + Ki J e 



ou 

- y* est la trajectoire de reference de la sortie, 

- e = y — y* est I'erreur de poursuite, 

- Kp, Kj sont les gains usuels. 



Kj = conduit au correcteur proportionnel intelligent, ou 
iP, 

F-y*+Kpe \ 

u = (3) 

a 

Remarque 1: Si = 1 en (IT|), on se ramene avec ([2]) ou 
(121) a la stabilisation d'un integrateur pur. D'ou le reglage 
facile des gains. 

III. MiSE en (EUVRE 

A. Estimation de F 
Reecrivons ^ sous la forme 

F = -au + y* - Kpe - Kj J e (4) 

On attenue les bruits entachant les mesures en integrant[f| 
les deux membres de ^ sur un court laps de temps S. II 
vient : 

-Flpprox = ^ J (^-au + y* - Kpe - Kj J dr 

ou ^approx est une approximation, constante par morceaux, 
de F. Get estimateur s'implante facilement sous la forme 
d'un filtre lineaire discret. 

Remarque 2: II y a quelques situations ou il peut etre 
interessant d'estimer a en ([T]). Voir, a ce sujet, ^36J, [41] . 

B. Trajectoire de reference pour la sortie 

On suppose le systeme d'entree u et de sortie y a de- 
phasage minimal|f|. On peut alors choisir une trajectoire de 
reference pour y satisfaisant les necessites du systeme. 

IV. Pourquoi l'ordre 1 suffit-il ? 

Prendre i' — 1 en. ^ suffit jusqu'a present dans toutes 
les applications^. Pourquoi cette heureuse propriete qui 
simplifie passablement I'implantation numerique du sans- 
modele ? On I'explique par les frottements. II leur corres- 
pond la presence de la derivee premiere y dans I'equation, 
inconnue, du systeme, qui evite I'apparition d'une boucle 
algebrique avec v = 1. lUustrons cette these avec le systeme 
lineaire du second ordre, a coefficients constants, 

X -\- cx + Ax — u (5) 

ou cx represente des frottements. Les figures (U [2] four- 
nissent des simulations satisfaisantes avec un regulateur 
iPI ; le choix des valeurs des parametres est c = 3, a = 1, 
Kp = 16, Ki = 25. Une degradation considerable des per- 
formances, comme on le voit sur la figure|3l se produit avec 
un oscillateur harmonique, c'est-a-dire avec c = en ([5]). 

3. Les bruits, consideres comme des fluctuations rapides, sont at- 
tenues par des filtres passe- bas, dont I'integrale est un exemple simple. 
Voir [16j pour une explication mathematique. 

4. L'ignorance des equations gouvernant le systeme interdit toute 
verification mathematique de la nature du dephasage. L'appreciation 
de cette nature repose done sur une bonne connaissance experimentale 
du comportement. Voir 1421 pour une premiere approche du dephasage 
non minimal dans le cadre du sans-modele. 

5. Sauf en [9], oix u = 0. 





Fig. 2 

COMMANDE DU REGULATEUR iPI. 



V. EfFICACITE DBS PI CLASSIQUES 
A. Discretisation 

Associons au correcteur PI classique 

u{t) — kpe{t) ~^ J e(T)(iT 

sa forme « vitesse » : 

u(t) = kpe{t) + kie{t) 
Une discretisation elementaire fournit : 



u{t)-u{t-h) _ fe{t)~e{t-h) 
h " I h 



(6) 



+ k,e{t) (7) 



ou h > est « petit » . 

Si = 1 en ([T]), remplagons F par y{t) — au(t — h) pour 
les besoins de I'implantation numerique. Alors, ([3]) devient 

= u{t -h)~ '^'^ "5' ' + — <t) (8) 
ha a 

(O et ([5]) deviennent identiques si, et seulement si, 



fc k--^ 
atl all 



(9) 



Remarque 3: Cette equivalence entre PI et iP n'est plus 
valable en temps continu, comme on le voit si /i — 0. 



B. Explication 

Les calculs precedents, qui s'etendent aisement aux PID 
(voir [4]), expliquent pourquoi, grace a nos regulateurs in- 
telligents, les correcteurs PI et PID classiques, echantillon- 
nes, donnent des resultats satisfaisants pour des systemes 
complexes si leur reglage est approprie. Le benefice, consi- 
derable, de nos iPI est, comme deja dit, un reglage simple. 

Remarque 4-' Voir [4] pour plus de details. Rappelons 
que cette comprehension, nouvelle semble-t-il, des PID a 
ete degagee lors de I'implantation du sans-modele a un ve- 
hicule electrique [3]. 

VI. MODELE RESTREINT 

Soit le systeme masse-ressort : 

my = -kiy - k^y"^ + F{y) - dy + u (10) 



ou 



seule la masse to = 0.5 est connue ; 
ki = 3, fca = 10; d = 5 sont supposes mal connus, et 
on utilise les estimees : fci = 2, fcs = 7, d = 2.5 ; 
un frottement discontinue 



-0.3 -0.4(y + 0.25)2 -dy si z)>0 
0.3 + 0.4(?) + 0.25)2 -rfy si y<0 



Son expression mathematique, indispensable pour la 

simulation, est inconnue du correcteur. 
On utilise la platitude ([22], [23], [34], [44]) du modele li- 
neaire restreint 



my = -kiy - k^y^ - dy 



(11) 



pour determiner une commande nominale en boucle ou- 
verte : 

u* = my* + kiy* + h{v*f + dy* 

Les simulations des figures U] et [S] correspondent a u* . Po- 
sons u = u* -\- Au et introduisons un iPI par rapport a 
Au pour stabiliser e = y — y* autour de 0. On aboutit aux 
excellentes simulations des figures [5] et [T] obtenues avec du 
bruit. 

Remarque 5: Le modele restreint (|lip serf uniquement, 
comme on vient de le voir, a calculer aisement une tra- 
jectoire de reference en boucle ouverte, tache aise s'il est 
plat. 



Fig. 4 

COMMANDE NOMINALE DU RESSORT. 



Fig. 6 

COMMANDE iPI DU RESSORT DANS LE CAS BRUITE. 




Fig. 5 

Sortie avec commande nominale 



Fig. 7 

La SORTIE DU RESSORT DANS LE CAS BRUITE ET COMMANDE iPI. 



VII. AUTRES ILLUSTRATIONS NUMERIQUES 

A. Regimes glissants 

La figure |S] represente uii vehicule militaire autonome, 
emprunte a p9|. Son interet est du au non respect de la 
« matching condition » , bien connue des tenants de la com- 
mande par regimes glissants (voir, par exemple, |12| . 
|14| . |29) . |40|). ce qui entraine une mise en place penible 
de cette commande en |29], alors que le sans-modele s'ap- 
plique sans difficulte aucune. Avec une dynamique lineaire 
du type 

x^Ax + Bu + BiUJ (12) 

ou to est une perturbation pouvant dependre du temps, 
rappelons que cette condition correspond a. Bi = (3B, ou /? 
est une constante. 

L' angle 9b doit prendre une valeur donnee avec une com- 
mande u agissant sur et une perturbation issue 9i. Les 
raideurs et les coefficients d'amortissement des ressorts- 
amortisseurs en les articulations 1 et 2, donnes par la 
matrice A, sont assez elevees. Les equations dynamiques 
sont de la forme ([12]) : x = [6'fc, 9„i, 9mb, Ob, 9u,9i]'^, w(i) = 

),Frnb X sign{9b-9p)Y , 

ou 

- 9m '■ vitesse angulaire du moteur, 

- 9}] : vitesse angulaire du bras intermediaire, 

- 9p : angle de tangage du vehicule pris comme pertur- 
bation. 




Fig. 8 

Vehicule militaire autonome. 



^mb 



%) + 9p . 



- N : rapport de reduction du moto-reducteur, 

- 9i : vitesse angulaire du troisieme bras 9i , 

- ujim ■ frottement au niveau du moteur, 

- ujii : frottement au niveau du troisieme bras 6i, 
~ Tarn ■ couple applique au moteur, 

- Tai : couple applique a la charge, 

- : frottement, 

- Wi„i(t) = fdSign{Tam - Jmdp), 

- UJii{t) = fdsign{Tai - JlOp), 

- Jm, Ji ■ inerties. 

On garde les valeurs numeriques de [29] pour A, B. Bi, C 
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u en guise de ([I}, un correcteur iPI, ou 
Kp = 160, Ki — 6400, fournit de bons resultats, d'apres 
les figuresl^etlTUl Notre commande, fort simple, est robuste 
par rapport aux bruits de mesure, non consideres en |29] . 
Ella ne necessite pas, contrairement a [29], 

- la resolution d'inegalites matricielles, 

- la mesure de 9^, O^b et 0^. 





Fig. 10 

Commande. 



vient 



ou 



x{t) = Ax{t) + Adx{t -t)+ Bu{t) + Biuj{t) (13) 

le retard r varie entre et 9 ms (I'amplitude est 
moindre en [29], ou il varie entre et 3 ms), 
la matrice Ad, tiree de [53], est 



Ad = 
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La meme commande qu'au § IVII-AI fournit les resultats 
tout a fait satisfaisants des figures [TTj et [T^ Elle est, ici 
aussi, bien plus simple que celle proposee en ||29j. 
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Fig. 11 

Angle Si, 



-0.2' ^ ^ ' 
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Fig. 9 
Angle 65. 

B. Dimension infinie 
B.l Retards 

Continuous I'exemple du § IVII- Al Ajoutons, comma [29] . 
un retard sur I'etat, du aux frottements. Alors, ([12]) de- 



Remarque 6: Les termes retardes en (jl3p se retrouvent 
en F dans ([ij. lis sont pris en compte sans autre forme 
de proces par notre strategic de commande. Nul besoin 
done d'identifier ces retards, tache delicate, surtout s'ils 
sont variables. 

Remarque 7: Rappelons (voir, par exemple, [5], [38]) que 
le reglage des PID traditionnels repose souvent sur une 
approximation du systeme a commander par un systeme 
a retard. Ici aussi, ces retards sont inutiles, comme il est 
demontre en jTT], [18], avec nos correcteurs intelligents 




Fig. 12 

COMMANDE. 



B.2 Systemes a parametres distribues 
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Fig. 13 

ECHANGEUR DE CHALEUR. 



B.2. a Un echangeur de chaleur. Pour I'echangeur de cha- 
leur de la figure 1131 nos simulations numeriques utilisent 
une discretisation spatiale des equations aux derivees par- 
tielles 



Mici 



M2C2 



dT 

'dt 
dS_ 
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miCi 



m2C2 



dT 

dx 
dS_ 

dx 



AU{T 
AU{T- 



S) 
S), 



ou 



- <x < L; 

- T{x, t) est la temperature du fluide 1, et T(0, t) la com- 
mando ; 

- S{x,t) est la temperature du fluide 2, et S{L,t) la 
sortie ; 

- U est la conduction. 

Les valeurs numeriques des parametres sont empruntees a 
un echangeur reel f7J. 

B.2.b iPI et PID. ^ devient : 



S{L,t) = F+l%T{Q,t) 



(14) 



On impose a S{L,t) une trajectoire de reference douce al- 
lant de 270° a 600° Kelvin, en utilisant une fonction poly- 
nomiale du temps de degre 6. On evite ainsi I'excitation des 
modes rapides. Quelques tatonnements ont suffi pour arri- 
ver a la valeur 18 en (|14|) . Les gains de I'iPI, deduit de (fT4|) . 



sont Kp — 10, Kj — 17.36. La methode de Ziegler-Nichols, 
bien connue et assez lourde, fournit les trois gains du PID 
classique, construit pour les besoins de la comparaison : 

r(0, t) = y* + 1.8 e + J e + 0.75 e 

ou y* est la trajectoire de reference et e = y — y* ■ Les 
deux commandes fournissent de bons resultats comme le 
montrent les figures [TH [151 [11] et[T71 Notons, neanmoins, la 
plus grande sensibilite au bruit du PID, due a la necessite 
de calculer une derivee. 




Fig. 14 

S(L,t) AVEC iPI. 




Fig. 15 

CoMMANDE iPI. 



Certains phenomenes complexes sont inevitables : 

- les vitesses et les proprietes physico-chimiques des 
fluides sont variables ; 

- le vieillissement s'accompagne de rouilles et de depots 
qui peuvent modifier la conduction U et retrecir les 
tuyaux. 

Toute commande de qualite se doit de les maitriser. Sup- 
posons, done, les tuyaux partiellement bouches, avec des 
baisses de 30% du debit et de 50% de la conduction. Sans 
nouveau calibrage des correcteurs precedents, 

- les figures [T5] et [TH] montrent les performances inchan- 
gees de I'iPI ; 

- celles du PID se degradent notablement d'apres les 
figures [20] et[2ll 
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Fig. 16 

S{L,t) AVEC PID. 
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Fig. 17 

COMMANDE PID. 



B.2.C Autres commandes. H y a une litterature impor- 
tante sur la commande des echangeurs de chaleur, qui obeit 
a des points de vue tres varies : voirH, par exemple, [7], 
[10], [HI, IS], [30], [43], [45], [51], [53|. Les equations aux 
derivees partielles, qui s'imposent dans toute modelisation 
un tant soit pen fine, y sont evidemment presentes, mais 
ne semblent pas avoir produit une regulation exploitable. 
Apparaissent aussi la commande predictive, la commande 
robuste, la linearisation par bouclage, la logique floue, . . .. 
Notre correcteur iPI, simple a mettre en oeuvre et aux 
bonnes performances, semble plus prometteur. 

Remarque 8: On trouve en [21] une interessante appli- 
cation du sans-modele au thermique, qui souffre d'oscilla- 
tions, dues, a notre avis, a une derivation numerique im- 
propre. On y remedie par les methodes du § IIII-A[ par 
exemple. 

VIII. Conclusion 

Des publications ulterieures developperont les considera- 
tions ci-dessous. EUes ont pour ambition de provoquer une 
discussion « epistemologique », aussi ouverte que possible, 
en mathematiques appliquees. 

A. Automatique 

Quel sens gardent les themes traditionnels de la recherche 
en automatique si les theses de la commande sans-modele 

6. La liste ci-dessous n'a aucune pretention d'exhaustivite. 
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Fig. 18 

S(L,t) AVEC iPI. 
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Fig. 19 
Commande iPL 



sont « vraies0 » ? Que Ton songe par exemple, d'un point 
de vue strictement academique, a 

- la plupart des questions structurelles en lineaire et en 
non-lineaire de dimension finie, puisqu'il n'y a plus de 
mo deles ; 

- I'identification parametrique, car sans modeles point 
de parametres ! ; 

- la commande robuste et la commande stochastique, 
puisque I'iPI prend en compte les termes inconnus et 
les perturbations ; 

- la commande optimale, deterministe on stochastique : 
le caractere toujours arbitraire et artificiel du choix du 
criterelfl est accentue par la meconnaissance du mo- 
dele. 

Aujourd'hui, I'exploration des limites du sans-modele 
semble une voie a privilegier, en particulier par rapport 
a la dimension infinie, c'est-a-dire par rapport aux retards 
et aux equations aux derivees partielles. Les succes obte- 
nus, notamment avec des systemes oii la modelisation par 
equations aux derivees partielles est primordiale, comme en 
[33) . avec les installations hydroelectriques, et en [1], avec la 
regulation du trafic autoroutier, rendent cette investigation 
urgente. 

7. On entend par « verite » la confirmation des succes pratiques. 

8. En ingenierie, contrairement a la physique « pure », ou les lois 
fondamentales derivent de principes variationnels, les criteres utilises 
pour I'optimisation sont « bricoles ». 
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Fig. 20 

S{L,t) AVEC PID. 
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Fig. 21 

COMMANDE PID. 



B. Modelisation mathematique 

Une large part de la recherche universitaire appliquee 
souligne, par-dela I'automatique, la necessite de I'ecriture 
de modeles mathematiques fiables, dont les equations dif- 
ferentielles, ou aux differences, sont I'armature. C'est trop 
souvent un leurre ! Y renoncer 

- n'implique pas I'abandon des mathematiques ; 

- signifie, comme indique ici, et en [19J pour la finance 
quantitative, la poursuite de concepts et, done, d'outils 
theoriques nouveaux. 
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